Extended navigability of small world networks: exact results and new insights 
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Navigability of networks, that is the ability to find any given destination vertex starting from 
any other vertex, is crucial to their usefulness. In 2000 Kleinberg showed that optimal navigability 
could be achieved in small-world networks provided that a special recipe was used to establish long 
range connections, and that a greedy algorithm, that ensures that the destination will be reached, 
is used. Here we provide an exact solution for the asymptotic behavior of such a greedy algorithm 
as a function of the system's parameters. Our solution enables us to show that the original claim 
that only a very special construction is optimal can be relaxed depending on further criteria, such 
as, for example, cost minimization, that must be satisfied. 

PACS numbers: 89.75.Fb, 05.40.-a 



By endowing nodes with both well-connected local 
neighborhoods and long-range shortcuts, that dramati- 
cally reduce the distances to any other node, transport 
on small-world (SW) networks is both locally and glob- 
ally efficient Such feature has made SW networks 
appealing for several fields, such as social, computer and 
life sciences 0, S S HE H H • Unfortunately, taking full 
advantage of the small node-to-node distances requires a 
global knowledge of the system, that is, in general, not 
accessible. It is thus important to devise decentralized 
algorithms that rely only on local information and that 
are able to find good, although sub-optimal, routes from 
source to destination 0] . The analysis of a prototype de- 
centralized algorithm showed that the precise recipe used 
to establish the long-range shortcuts affects the ability of 
decentralized algorithms to navigate the networks [lfj. 
Here we provide an exact solution in any dimension for 
such problem, and, in light of recent insights into the 
properties of SW networks ll| and in contrast to the 
findings in we show that there is a broad range of 
SW networks that are optimally navigable. 

SW networks can be obtained from regular lattices in 
d dimensions by adding to every node, with probabil- 
ity q, a long-range connection to another node taken at 
random over the whole lattice |12j. The key feature of 
SW networks is that the average distance between any 
two nodes grows at most logarithmically with the linear 
size L of the lattice [lij]. In order to analyze how dif- 
ferent shortcut addition schemes affect t he p erformances 
of a decentralized algorithm, Kleinberg [lfj specified an 
additional rule for the selection of the long-range part- 
nerships: the probability s(l) that the shortcut added to 
a node ends at a node at euclidean (or lattice) distance 
I is a decaying power-law, s(l) = Nl~ a , where AT is the 
normalization over the whole lattice. The two key ingre- 
dients for the construction of the SW network are thus 
the shortcut addition probability q and the exponent a. 

The simple decentralized algorithm considered in [l~oj ] 
is of greedy nature: starting from a given node, at ev- 
ery step the algorithm chooses the edge with the end- 



node which is closer, in euclidean or lattice distance, to 
the selected destination; such scheme guarantees that the 
destination is always certainly reached. Using arguments 
from probability theory, it was possible to find, in d = 2, 
a lower bound for the average number of steps r(L) that 
are necessary to connect nodes separated by a distance 
proportional to the lattice linear size L: t(L) < L 13 , with 
(3 = (2 - a)/3 if < a < 2 and /3 = {a - 2)/(a - 1) 
if a > 2 (see Fig[TJ dashed lines). In the case a = 2 
the upper bound t(L) < (InL) 2 was found instead. As 
a consequence, the best choice for an optimized naviga- 
bility of a SW network using a decentralized algorithm 
would be a = 2 in d = 2 and, more generally, a — dm. d 
dimensions 1(J 13 1. 

In what follows, we derive the exact asymptotic be- 
havior of t(L) in any dimensions. We resort to the same 
implicit assumption already used in the probabilistic ap- 
proach in fl3| : we treat the algorithm as a stochastic 
Markov process, by looking simultaneously at all pos- 
sible network realizations [lj|. If at a given stage the 
message is at site i, at lattice distance di from the tar- 
get, at the next step it will surely be at a site j, with 
distance dj < di, due to the greediness of the algorithm. 
More in detail, if j is a nearest neighbor of i, then their 
connecting edge is chosen only if there are no shortcuts 
from i to sites k with dk < dj. If, instead, such a useful 
shortcut exists, the greedy algorithm chooses it over the 
nearest neighbor connection. In general, we can write 
the following recursive relation for the average number 
of steps from a site i to the destination 



(1) 



where the probability Pi-^j depends on the presence and 
greedy-usefulness of shortcuts. In equation (Q} we as- 
sume, following [l(J, that the time it takes to travel a 
network edge, be it a shortcut or a link of the underlying 
lattice, is 1. 

Equation ([I]) can be easily solved numerically, on a 
lattice, by recursion. However, it is more instructive to 
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take its continuous space limit, where a lattice site i is 
mapped onto a position r and the lattice spacing van- 

I 



ishes. Equation (p} in d dimensions then becomes (for a 
complete derivation see the auxiliary material (lEj ) 
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with e a short-lengthscale cutoff that avoids divergences 
when / — > (e was clearly not necessary on a lattice). In 
equation ©, J dil is the integral over the ci-dimensional 
hypersolid angle and it is parametrized, among others, 
by the azimuthal angle which covers only half of the 
hypersphere [l5j : q is the linear probability density of 
shortcuts. Clearly it takes no time to reach the desti- 
nation starting from itself, and consequently r(0) = 0. 
Given the complete isotropy of the problem, r(r) depends 
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only on r = \r\. 

By setting the scaling function form 

r(r) = K^fiKr) 



(3) 



(the value of K to be discussed below case by case), it 
is possible to obtain, after some manipulations 15|. an 
integral equation for the derivative of f{x), 



y — cos ( 



^Jy 1 + 1 - 2ycos6»" 



/' (xy/y 2 + l-2ycos0) dy 
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where I(y,Ke/x) = /* {Ke/x + z)~ a z^dz. To ob- 
tain the asymptotic behavior of f(x) we rely on sim- 
ple considerations. Since r(r) cannot decrease with the 
distance r from the target, f{x) is non-decreasing and 
therefore f'(x) is non-negative. Moreover, T(y, Ke/x) is 
positive and thus from ^ we obtain f'(x) < 1. It is 
then possible to show that if a < d + 1 , asymptotically 



f'(x) depends only on a and f'(x) 
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stead a > d + 1, f'(x) is not universal anymore [15| and 
f'(x) — > 1/[1 + c(a, e, d)g], where c(a) is a constant de- 
pending only on d and a. 

In the following we analyze the five cases < a < d, 
a = d, d < a < d+l, a = d+1 and a > d + 1 separately: 
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d: In this case the normal- 
= (d - a)/L d - a and X = 
a) £-(<*-«)/(<*+!-«) to the i ea ding 

Since /'(a;) - l/a;( d+1 - Q ), /(x) 



converges asymptotically to a constant and con- 
sequently t(L) - L (d-a)/(<i+i- Q )/ g i/(d+i- Q ) for 

large values of L. 

a = d: The normalization constant is M = In L 
to the leading order in L. This implies f'(x) ~ 
1/x, so that r(r) ~ (InL) /g for large values of L. 
This results coincides with the upper bound found 
in 0. 

d < a < d + 1 : The normalization does not de- 



pend, to the leading order, on L. In this case 
f'(x) ~ is non-integrable and there- 



fore f{x) diverges as 
r(L) - L Q - d /g. 
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and asymptotically 



• a = d+l: Again the normalization does not depend 
on L. After some calculations [15| it is possible to 
show that r(r) ~ L/{qhxL). 

• a > d + 1: Once more the normalization does not 
depend on L. It is easy to show that f'(x) — > 
l/[l+c(a, e, d)g] and thus r(i) ~ i/(l + c(a, e, d)g) 
asymptotically. 

The above results are summarized in TablelH and one ma- 
jor conclusion that can be drawn at this stage is that a 
greedy algorithm can significantly outperform the simple 
lattice distance only if the shortcut length distribution 
has diverging first moments in the L — > oo limit. The 
analytical predictions, shown in Fig[T] (red solid lines) 
for d = 1,2, are compatible with the lower bounds found 
in [nj]. We have also numerically solved equation (JlJ on 
a lattice, and using lattice distances, in d = 1,2 and the 
estimates for the exponent (3 agree with the analytical 
predictions (Fig[T] blue stars), apart from small discrep- 
ancies due to the not-yet fully achieved asymptotic limit, 
an effect that is expectedly more important in higher di- 
mensions. Furthermore, in d = 1 we have verified that, 
as long as a < 2, the r(r) curves for different values 
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of q and of L, but for the same value of a, do collapse 
onto each other asymptotically once Kr(r) is plotted as 
a function of Kr with the appropriate values of K (see 
Figf2]). This result confirms the asymptotic universality 
of f'(x) when a < d + 1. 

The results described above are necessary but not yet 
sufficient to thoroughly address the navigability of SW 
networks. Indeed, changing the value of a while keeping 
the value of q fixed, as in [loj], is only one of the pos- 
sible ways to compare the navigability of SW networks. 
As a matter of fact, it has been shown that, as long as 
< a < 2d, it is possible to make a <i-dimensional lat- 
tice small by letting the shortcut probability q depend 
on a [111 ]. In particular, the shortcut probability q mark- 
ing the crossover from the euclidean to the small world 
regime has the form q ~ L~ d if < a < d, q ~ In LL~ d 
if a = d and q ~ L a ~ 2d if d < a < 2d. Thus, keeping 
q fixed for different values of a would pit against each 
other networks that are intrinsically differently small. 

To alleviate such bias, a more appropriate comparison 
should take into account how q must change with a. In 
order to do so, we can start from a given value of the 
shortcut probability q{d) at a = d, where purportedly 
navigation is easier. Then, as a moves away from d, 
we let this value change according to q(a) = q(d)/\nL 
if a < d and q(a) = q{d)L a - d / lnL if d < a < 2d. 
Since our derivation of the asymptotic behavior of t(1) 
also provides the precise form of the g-dependent prefac- 
tors, we can obtain the correct expected asymptotics of 
t(L) when q is allowed to appropriately change: t(L) ~ 

(L d - a lnL) 1/{d+1 ~ a) if < a < d, r(L) ~ (lnL) 2 if 
a = d, t(L) ~ lnL if d < a < d + 1, r(L) ~ const if 
a = d+l and r(L) - L-[<*-(d+i)} if d+1 < a < 2a >. The 

last result is apparently paradoxical, because it predicts 
the navigation time to become smaller for larger systems. 
Actually, using the results obtained in [llj, it is possible 
to show that if q(a) = q(d)L a ~ d / lnL, the average path 
length of a small- world actually decreases with L for fixed 
a, even faster than using a greedy algorithm [ljj ]. It is 
thus not surprising that also within the greedy framework 
larger values of a, with the corresponding increase of q, 
lead to shorter navigation times. Considering all the re- 
sults summerized in Table U we conclude that if we take 
the necessary change of q into account any d < a < 2d 
outperforms the case a = d. 

To gain a definitive insight into the properties of such 
class of SW networks, we finally consider the total short- 
cut length per node C = q(a, L)l(a, L) (l(a, L) being the 
average shortcut length), for different values of a. C can 
be associated to the amount of resources (cost) needed 
to set up such networks. We have C = q(d)L/h\L for 
0<a<d+land£ = q(d)lL a ~ d / lnL for a> d+l. 
We observe that the value of L is the same for all 
< a < d + 1, and lower than that for d + 1 < a < 2d. 

The results presented above reveal a new picture for 



the navigability of SW networks, which is richer and more 
variegate than previously outlined: navigability depends 
on several parameters, and the optimal choice depends 
on the criteria that have to be satisfied. If the average 
number q of shortcuts per node must be kept constant, 
at the cost of making networks characterized by differ- 
ent values of a differently small, then our exact results 
confirm, while making more precise, the findings of [loj ]. 
If on the other hand a fairer comparison of small-world 
networks is desired, then the necessary change of q as a 
function of a must be taken into account and any a > d 
outperforms a — d, with optimal navigability reached 
for a = 2c?. If finally the amount of resources needed to 
set up the networks becomes a crucial factor, the whole 
range a £ (d,d+l) is where navigability is optimal and 
cheapest. 

Outlining the dependence of navigability on various 
criteria, as highlighted in this work, is important in or- 
der to understand the architecture of networks where the 
small-world feature is believed to be crucial. Indeed, it 
can be expected that virtual connections, such as hyper- 
links between web-pages, will bear almost no signature 
of cost-effects, but might be limited in number to a few 
units; contrary- wise, the number of connections per neu- 
ron in the brain can be extremely high, but establishing 
and maintaining them is surely resource expensive. Thus 
different systems might have achieved different layouts 
to enhance their navigability while respecting given sets 
of constraints. Keeping in mind all these criteria can be 
of paramount relevance given the growing interest in the 
greedy navigability of real networks [lij ]. 

The authors would like to thank the Swiss National 
Research Foundation (200020116286) for financial sup- 
port. 
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r(L,q) 


r(L,q(L,a)) 


< a < d 


d-a 
Ld+l-a 


( L d ~ a \ d+l-oc 
\ l»i ) 


a — d 


(lni) 2 


(InL) 2 


d < a < d+ 1 


L a ~ d 


InL 


a = d+ 1 


L/lnL 


const 


d + 1 < a < 2d 


L 


L -[ a -(d+i)] 



TABLE I: Asymptotic navigation times as a function 
of the exponent a. The column t(L, q) reports the navi- 
gation time at constant shortcut density, whereas the column 
r(L,q(L,a)) reports the navigation time when the shortcut 
density is allowed to change so to keep the system in the 
small-world regime. 
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FIG. 1: Asymptotic exponents of the greedy distance. 

Comparison between the lower bounds found in the original 
reference (black dashed lines) and the exact solution presented 
here (red solid lines) in (a) d = 1 and (b) d = 2. The exact 
solutions are verified by numerical results (blue stars). The 
shadowed region (a > 2 in d — 1 and a > 4 in d — 2) is 
outside of the small-world regime. 




Kr 



FIG. 2: Universality of greedy distances. The numerical 
results confirm the ansatz r(r) = K -1 f(Kr) in d = 1. The 
collapses obtained once Kr(r) is plotted as a function of Kr 
for a fixed a but for different (g, L) pairs show that asymp- 
totically (Kr — > oo) f(Kr) is a universal function that only 
depends on a. 

Supplementary Material 

STOCHASTIC PROCESS FORMULATION 

Kleinberg's problem can be formulated as a stochastic process when analyzed on the whole small-world network 
ensamble. We consider a walker that starts from site i at distance di from the destination. At the next step it will 
move differently on every network realization in the ensemble. Thus it will be on a site j at distance dj from the 
destination with a probability q, that is the probability that site i has a shortcut, times the probability that there is 
a shortcut from i to j, times the probability such shortcut can be used by the greedy algorithm, i.e. that dj < di. 
When j is a nearest neighbor of i (with dj = di — l), the probability of the walker to be in j is given by the probability 
that i does not have any shortcuts or that it does but that they end at distances from the destination larger than dj . 
In that case the algorithm chooses an edge of the underlying lattice. It is then possible to formulate the evolution of 
the probability distribution of the walker's position as 

P k (t + l)=J2p^kP J (t) (5) 

j 

with initial conditions -Pfc(O) = (5^ if the walker starts from i. Given the nature of the greedy algorithm, the destination 
will be reached with certainity in a number of steps that is at most di. 

Once the stochastic formulation is established, it is then easy to derive that the average trip time Ti from site i to 
the destination can be expressed recursively using the average trip time from every other node closer to the destination 
than i: 

T l =J2'P^j( T J+ 1 ) ( 6 ) 

j 

where the 1 in the parenthesis indicate the time it takes to go from site i to site j if there is a usable edge (shortcut 
or otherwise) from i to j, which, as explained above, occurs with probability Pi^j. Equation ^ can be easily solved 
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with the aid of the computer, but it is extremely demanding on time and memory resources for lattices with dimension 
d > 2 if the asymptotic limit has to be reached. 



CONTINUOUS SPACE FORMULATION OF EQUATION © 

The continuous space formulation of equation ([6]) is easily obtained as follows. For reasons that will be clear in the 
following, shortcuts do not connect lattice sites, but rather the shortcut associated to a site go from the lattice edges 
connecting it to its neighbors to a lattice node. A small-world network obtained in this way is clearly as small-world 
as a network obtained with shortcuts from site to site. 

Without affecting the generality of the results, from now on we shall assume that the destination is the origin. 
We also assume spatial isotropy of the process. Say that the walker starts from a point r, at distance r from the 
destination. Since it obeys a greedy scheme, it will always try to approach the destination as much as possible. Thus, 
in the absence of any useful shortcut, it will move along the ray defined by r. Let's examine what may happen while 
it moves along the ray for an infinitesimal distance dr. The point at distance r — dr will be reached only if 1) no 
shortcuts are encountered, with a probability 1 — qdr (where q is now the linear density of shortcut probability, which 
also explains why the starting points of shortcuts are more appropriately associated to lattice edges in this context), 
or 2) if shortcuts are encountered (probability qdr) but none of them are useful. The expected time to reach the 
destination, in d dimensions, can then be written as 



r(r) 



(1 - qdr) + qdr ^1 - J dQ J 



JV(e + j)~ a r~ 1 dz 



t{t - dr) + — ) + 



(7) 



where the dots indicate the contribution of useful shortcuts encountered along dr, that we will discuss below. J\f(e+l)~ a 
is the shortcut length probability distribution and e is a short-distance cutoff introduces in order to avoid divergences 
when I — ► if a > d. 

The double integral in equation ([7]) is the probability for a shortcut encountered between r and r ~ dr to fall inside 
a sphere of radius r, thus closer to the destination than r. It is composed of an angular integral over all the angles 
describing a (hyper)spherical parametrization of the space centered on r (with 9 S [0, 7r/2] in the integral over dQ), 
and of a radial part I centered on r. The triangle defined by r, I and r — I identifies a two-dimensional plane in any 
dimension. Given the isotropy of the space, there is full rotational invariance around the ray r, hence 9 is the only 
angular variables of the integral over dfi that affects the radial integral. The reason for this peculiar choice of the 
parametrization will be explained below. 

The term t(t — dr) + — is the expected time to reach the destination from a distance r — dr augmented of the time 
to reach r — dr from r, which we assume to be proportional to dr and to depend on an intrinsic velocity v. In what 
follows, we fix v = 1. 

Equation ([7]) is not complete yet. The term representing the contribution from useful shortcuts is given by: 



-(r) 



qdr / dfl 



2r cos ( 



J\f(e + l^V 1 - 1 (t (y/l 2 + r 2 -2/rcos0) + ^\ di 



(8) 



In complete analogy with our assumption, in equation ([6|) , that it takes exactly the same unitary time to travel along 
a lattice edge or a shortcut, the time to reach r — dr from r is the same in (0 and (|8|) (with v — 1 in the latter). 

By dividing both terms by dr and taking the limit dr — > 0, we obtain the continuos space limit presented in the 
main text: 



r'(r) = l-gTVr(r) / dQ 



2r cos 8 



(e + iyr^di + qjV / dn 



2r cos 8 



{e + l)- a l 



a jd—l„ 



l 2 + r 2 



2lrcos9) dl 



(9) 



FIG. 3: Graphical representation of the relations between r, r' and I. 



RESCALING OF EQUATION © 

When we introduce the assumption r(r) = K^ 1 f(Kr) the left hand side of equation © becomes r'(r) = f'(Kr). 
After some manipulations, we can write: 

/> />2rcos# 

f'(Kr) = 1- qNK-^ d+1 -^f{Kr) / dVl / (Ke + Kl)- a (Kl) d - 1 d(Kl) + 



+ q AfK- (d+1 -^ J dn J (Ke + Kl)- a f (yKH 2 + K 2 r 2 - 2(Kl)(Kr) cos 0) (Klf^diKl) .(10) 
A simple change of variable y = Kl/ Kr in the integral leads to 

/p2 cos 9 / -rr \ —ot 

dn J a \K~r +V ) yd ~ ldy + 

p p2 cos B / jy- \ —a 

+ qJVK- (d+1 - a) (Kr) d ~ a J dnj \j^. +v ) / [Kry/y* + 1 - 2y cos 6>J y d ~ l dy . (11) 
Now we can set the value of K according to qMK~^ d+1 ~°^ — 1. Once the corresponding expression of K is introduced, 
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equation (llip depends only on Kr. We can, then, change the general variable to x = Kr and obtain the equation 



/•2 cos 6 / Ts- 

f( x ) = i-x d ~ a f(x) I dnj^ (— + y) v d ~^y 



x d - a J dfl J \— +y ) f (Wj/ 2 + 1-2j/cos0J if- 1 Ay . (12) 



TRANSFORMING THE INTEGRO-DIFFERENTIAL EQUATION INTO AN INTEGRAL EQUATION 

By calling I(y,Ke/x) the primitive of (— + y) y^ 1 , equation (fTSj) becomes 
f\x) = l-x d - a f{x) JdSl [1(2 cos 6, Ke/x) -1(0, Ke/x)] + 

p <>2cos0 f Ke \~ a / \ 

+ x d - a Jdnj (~ + y ) y d ^f{ x Vy 2 + i-2y cos 9) d y . (13) 

Integration by parts the last integral in equation (jT3J) gives 

f'{x) = l-x d - a f(x) [ dO [7(2 cos 6, Ke/x) -1(0, K e/x)} + 



x d ~ a I dSl{I(2cos0,Ke/x) - 1(0, Ke/x)] f(x) 

dSl I( y ,Ke/x) —- y = f (x^/y2 + l-2yco S e)dy (14) 

J Jo \/y 2 + 1 - 2y cos 6 v > 



which simplifies to 

^2 cose 



f'{x) = l-x d+1 ~ a [dQ [ "' I(y,Ke/x) — y C ° Sd f'(xy/y* + 1 - 2y cos ^) dy . (15) 
J Jo V?/ 2 + 1 - 2y cos 6 v > 

This is an integral equation for fix). The properties of the kernel I(y,Ke/x) are simple: it is always positive 
because it is the integral of a positive function, and 1(0, Ke/x) ~ (Ke/x) d ~ a . Thus, in the asymptotic limit x — > oo, 
I (y, Ke/x) develops an integrable singularity for a < d+ 1 and it converges to a well-behaved kernel. In what follows 
we can safely disregard the short-distance cutoff. For a > d + 1, the singularity is non-integrable and it will require 
a special treatment. 

Equation (fT5|) shows that, since we can safely take the Ke/x — > limit when a < d + 1, f'(x) must obey an 
asymptotically universal form that depends only on a. We have numerically verified that this is indeed the case (see 
main text). In the case a > d + 1, instead, the Ke/x term dominates and f'(x) will not be universal any more. 



ASYMPTOTIC BEHAVIOR OF f (x) 

We want to determine the asymptotic behavior of f'(x). Clearly r(r) must be a non decreasing function of r, since 
it obviously takes longer to reach the destination starting from increasingly farther distances. Thus, f'(x) is non 
negative. Consequently, equation (fT5)) tells us that f'(x) < 1. We can, then, rewrite equation (fT5|) as 



fix) = 1- / dfi 



2 cos 6 



I(y,Ke/x) 



y — cos ( 



1 — 2y cos 



d+2-a 



1 — 2y cos i 



d+l-a 



1 — 2y cos ( 



(16) 
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The case a < d + 1 

If a < d+1 it is easy to show that the only ansatz for the behavior of x d+1 ~ a f'(x) that does not lead to contradicting 
conclusions is that 



"fix) 




I(y,Ke/x) 



y — cos ( 



y 2 + 1 — 2y cos 9 



(17) 



so that f'(x) ~ C/x d+1 ~ a . Indeed, should f'(x) decrease faster than that, then the whole integral would become 
increasingly negligible and f'(x) would be asymptotically equal to 1, against its fast decay assumption. Should 
f'(x) decay slower than C/x d+1 ~ a , then the integral would asymptotically diverge and f'(x) would become negative, 
implying a non-physical decrease of r(r) for larger distances. 



The case a — d + 1 



In this case K is not set by any equation, hence we can choose K = 1. This implies that in the rescaling of the 
continuos limit we can simply set r(r) = f(x) and r — x. The integral equation becomes (already in the limit when 
x — > oo) 

/> />2 cos 6 t\ 

f'(x) = 1-qMUn / (%) ./" = = /' (x^y 2 + l-2ycos6)dy . (18) 

J Jo ^Jy 2 + l- 2ycos6 v > 

Since asymptotically J 2cos I(y, ajx)dy ~ lnx/e, we have f(x) = 1/[1 + qc(a, d) Lax]. Thus r(r) ~ g7 v e ( a d) 



The case a > d + 1 



If a > d + 1, the kernel diverges as (x/Ke) a d for a; — ► oo. We thus divide and multiply it for its integral between 
and 2 cos 6* to obtain 



f(x) = l-x 



d+l-a 



2 cosfl 



I(y 7 Ke/x)dy 



y — cos ( 



2cosfl 

<\il I l(y,Ke/x)- 

o v 2/ + 1 — 2y cos y 



/' (Vt/ 2 + l-22/cos0j dy 
(19) 



where X(y, Ke/x) — I(y, Ke/x)/ J Q 2 cot ' e I(y, Ke/x)dy . Interestingly, T(y, Ke/x) — > <5(y) as a; — > oo and thus equa- 
tion ([±9]) reduces to 



/'(as) = 1 - a 



.d+l-c 



2 cos 6 



I(y,Ke/x)dy 



fix) 



(20) 



Since J Q cos I(y, Ke/x)dy ~ (a;/i4re) Q_rf_1 , from equation (|20|) we obtain that /'(a;) ~ 1/(1 + c(a, e, d)q), where 
c(a, e, d) is a constant that depends only on a, e and d. Thus, /(a;) ~ a;/ (1 +c(a, e, and r(r) ~ r/ (1 + c(a, e, c?)q). 

This result can be derived on a lattice by simple arguments. When a > d+1, the average shortcut length I is 
finite, so that shortcuts simply redefine the elementary length-scale of the system. The typical step length is, then, 
hyp = 1 ' {(1 ~ 0) + ?[1 ~ c(a, e, d)]} + Z • gc(a, e, d) = 1 + (/ — l)gc(a, e, d), where c(a, e, d) is the probability that, if 
a shortcut is present, it can be used by the greedy algorithm. I > 1 because on a lattice there are no shortcuts of 
unitary length. The average time to reach the destination, given that it increases linearly, will be r/lt yp . 



GREEDY AND GLOBAL PATH LENGTH IN THE CASE a> d+1 AND q(a) = q(d)L a ~ d / In L 



As outlined and numerically verified in [Petermann and De Los Rios, Phys. Rev. E 73, 026114 (2006)], the average 
path length < L > for small world networks of linear size L built according to the present recipe (given a and q) 
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obeys the scaling relation 

<r L ->= L*T. ( ^ 

where L* = q-^ d if (a < d) and L* = g-V(2<*-°0 if (d < a < 2d). The function T a {x) is 



<L>=L*tJ±-\ ,21) 



= n vr fi! (22) 

I (lnx) 7 if x 3> 1 

where the x » 1 asymptotic behavior is a power of a logarithm (see the above mentioned reference and also [Kosmidis, 
Havlin and Bunde, Europhys. Lett. 82, 48005 (2008)]). 

Using q — q(d)L a ~ d / InL (for a > d) in equation (f2Tj) it is possible to obtain 

< L >= L-^{\nL)^T a ( L2d ~'* 1 | . (23) 

V(lnL)^ J 

Thus, since T a increases at most logarithmically, the average path length decreases asymptotically as a power of L 
as a increases. The logarithmic decrease of the greedy path length for a > d+ 1 is thus compatible with this behavior 
and, expectedly, the greedy path length remains larger than the global one. 



